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ABSTRACT 

The Stephan's Quintet (SQ, HCG92) was observed with the Far-Infrared Sur- 
veyor (FIS) aboard AKARI in four far-infrared (IR) bands at 65, 90, 140, and 
160 fim. The AKARI four-band images of the SQ show far-IR emission in the 
intergalactic medium (IGM) of the SQ. In particular, the 160 //m band image 
shows single peak emission in addition to the structure extending in the North- 
South direction along the shock ridge as seen in the 140 /zm band, H 2 emission 
and X-ray emission. Whereas most of the far-IR emission in the shocked region 
comes from the cold dust component, shock-powered [C II]158/xm emission can 
significantly contribute to the emission in the 160 /j,m band that shows a single 
peak at the shocked region. In the shocked region, the observed gas-to-dust mass 
ratio is in agreement with the Galactic one. The color temperature of the cold 
dust component (~20 K) is lower than that in surrounding galaxies (~30 K). We 
discuss a possible origin of the intergalactic dust emission. 

Subject headings: ISM: structure — galaxies: clusters: individual (Stephan's 
Quintet) — galaxies: interactions — infrared: ISM 
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Introduction 



The track of the galactic scale transmigration of gas and dust is indispensable to un- 
derstand the gala xy evolution. S tephan's Quintet (SQ, HCG92) is the well studied compact 
group of galaxies ( IHickson Ill982l ) with the disturbed intergalactic medium (IGM) as shown in 
Fig. [TJ The intruder galaxy NGC 7318b is currently colliding with the IGM at a relative veloc- 
ity of ~ 1000 km s _1 and causes a large-sc ale shock front and IGM starb ursts (SQ-A and SQ- 
B). The shock front was first discovered by lAUen and Hartsuiker I (119721 ) with radio, and sub 
sequently detected in the X-ray em i ssion from shock-heated ( ~ 6 x 10 6 K) gas (Pietsch et al. 



1997; Trinchieri et al. 2003 



(2006) found 



20051 : lO'Sullivan et al. Il2009h . lAppleton et al. 
powerful H 2 rotational line emission from warm (~ 10 2 -10 3 K) molecular gas in the center 
of the shock ridge. The significant discovery is the extremely large equivalent width (EW) 
of the H2 line emission. This may indicate that far-infrared (IR) fine-structure lines such 
as [C II] 158 fim line also show extremely large EW by the shock and have a significant 
contribution to far-I R luminosity. However, the fact has yet to be revealed observationally. 
Cluver et al. I (120101 ) made H 2 maps of the SQ and found that the spatial distribution of 



the H2 emission was similar to those of radio an d X-ray emissio n in th e shocked region. To 
explain the co-existence of both hot and H 2 gas, iGuillard et al. I ( 120091 ) proposed a model of 
the shock in the inh qmogeneous gas medium. While H 2 molecules are likely to be produced 
on the grain surface (IGould fc Salpeter Ill963l). direct association of the H 2 gas with dust has 



yet to be found in the shock ridge. Recently, iNatale et al. I ( 1201 01 ) reported dust emission in 
the SQ with Spitzer observations. While the map at 160 /xm shows the significant presence 
of cold dust in the shock ridge, it is difficult to make direct comparison of the diffuse 160 /im 
map with the H 2 emission map because of the effect of bright emission from the surrounding 
galaxies. 



The Far-Infrared Surveyor (FIS) (IKawada et al. 1120071 ) aboard AKARI ([Murakami et al. 

20071 ) has four far-IR bands with the central wavelengths of 65, 90, 140, and 160 /xm, 
and achieves high sensitivity with relatively high spatial resolution. Finer allocation of 
AKARI /FIS 4 bands can provide the spectral information to constrain both the dust tem- 
perature and the [C II] 158 //m line emission and is the distinctive advantage compared to 
Spitzer /MIPS observations. In this paper we report the observation of the SQ with the FIS, 
and discuss the origin of the far-IR emission in the IGM. 



2. Observations and Data reduction 



The SQ was observed i n part of the AKARI mission program "ISM in our Galaxy and 
Nearby Galaxies" (ISMGN; iKaneda et al~ll2009f ) on 2007 Jun 18 using the FIS01 observa- 
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tion mode (Observation ID: 1402238-001 ). Details of the F IS inst rument and its i n-orbi t 



performance/calibration are described in lKawada et al. I (120071 ) and lShirahata et al. I ( 120091 ). 



The FIS was operated in the photometry mode with the four bands: N60 (65 /mi), WIDE-S 
(90 /mi), WIDE-L (140 /mi), and N160 (160 /mi). The four-band data were simultaneously 
obtained in a pointing observation for an area of about 10' x 20' around the SQ. The FIS 
data were processed with the AKARI official pipeline modules (version 20070914). Finally, 
four-band images were created with grid sizes of 12.5" for the WIDE-L and N160 bands and 
7.5" for the WIDE-S and N60 bands. The widths (FWHM) of the point spread functions 
(PSFs) are ~ 45" for the WI DE-L and N160 bands and - 30" for the WIDE-S and N60 



bands (IShirahata et al. II2009I ). The four-band images are presented in Fig. [2j The images 
are smoothed with boxcar kernels with a width of 25" for the WIDE-L and N160 bands, and 
15" for the WIDE-S and N60 bands. 

To derive the spatial variation in color temperature, the spatial resolutions of the original 
WIDE-S and N60 images are reduced to match those of the WIDE-L and N160 images by 
convolving the former images with a Gaussian kernel. The images are then resized with the 
common spatial scale among the four bands, 12.5"/pixel. The four-band flux densities at 
each image bin are derived after subtracting the sky background. The background levels were 
estimated and subtracted using nearby regions of the blank sky, which were observed in the 
beginning and at the end of the scan. Color corrections were also applied for the obtained 
flux densities by assuming a modified blackbody spectrum with the emissivity power-law 
index of unity and the temperatures of 20 K and 30 K for the shocked region and galaxies, 
respectively. Aperture corrections were not applied to the flux densities because we only 
consider diffuse emission components. At the centers of the shocked region, NGC 7319, and 
NGC 7320, we derive the flux densities by integrating the surface brightness within circular 
apertures. Figure [3^a) shows the resulting spectral energy distributions (SEDs) at these 
positions. Note that the fluxes in Fig. [3£a) are not the fluxes for the whole areas of the 
shocked region and the galaxies. The systematic errors for the flux densities are estimated 
to be -10 % for the WIDE-S and N60 bands, ~20 % for the WIDE-L band, and -25 % for 
N160 band, which include the uncertainty due to radiation effects. The relative uncertainties 
for the flux densities are estimated to be below 5 % for the four bands. 



3. Results 

In Fig. El N60 and WIDE-S images show emission from NGC 7319 and NGC 7320, 
and from the SQ-A and SQ-B (Fig. 121(c)). Although the WIDE-L image shows a structure 
extending in the North-South direction along the shock ridge, the N160 image clearly shows 
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single peak emission in addition to the structure, but does not indicate any features associ- 
ated with the individual galaxies. In the shocked region, the background-subtracted surface 
brightness in the N160 band is larger than 9 MJy sr^ 1 , which is comparable to 5cr sky back- 
ground noise. For the first time, AKARI observations reveal the dramatical change in the 
spatial distribution of far-IR emission at wavelengths ar ound 140-160 urn. The difference in 
images between the N160 band and the MIPS160 band (INatale et al. 1120101 ) may come from 
the fact that the spectral response of the N160 band covers longer wavelengths than that of 
the MIPS160 band; wavelength coverage at 10 % of the peak responsivity is 140-196 fim for 
the N160 band (AKARI FIS Data Users Manual ver.1.3) and 129-184 /im for the MIPS160 
band (MIPS Instrument Handbook). Figure [2]^f) shows that the spatial distribution of 160 
jum emission is quite similar to that of X-ray emission in the shocked region. In fact, X-ray, 
radio, and H 2 emission in the shocked region all show structures associated with the shock 
ridge and the bridge extending eastward from the shock ridge to NGC 7319. Our result 
clearly shows the spatial correlation of the far-IR emission with these structures. 



3.1. Difference in images between N160 and WIDE-L bands 

The cold dust temperatures in the shock and surrounding galaxies are about 20-30 
K. Because the flux ratio between the N160 and WIDE-L bands is not very sensitive to 
the dust temperature, the difference in the spatial distribution of far-IR emission between 
the two bands is hard to explain only by the presence of the cold dust component. An 
alternative possibility is a contribution from the shock-powered [C II]158/xm line emission 
to t he N160 band. For the SQ, the wavelength of [C II]158/xm line is redshifted to 161 



jum (IHickson 1119921 ). Thus, its contribution to the N160 band is larger than to the WIDE-L 
band. To estimate the contribution from the [C II] line, the WIDE-L flux density Fwl is 
simply subtracted from the N160 flux density -Fni60 by assuming that the dust emission 
is constant in units of F„ over the two bands (Fig. E£b)). By using the same aperture as 
that in Fig. [2^f), the subtracted flux density _F su b is estimated to be 40^22 mJy. The error 
includes the uncertainty that F WL to be subtracted must be increased by 38 %, if the SED of 
the dust emission is considered as the modified blackbody spectrum whose temperature and 
amplitude are determined from the data at 90 and 140 /im. The uncertainty is comparable 
to the subtracted sky background noise. The [C II] luminosity surface density £l [cii] [erg 
sec -1 kpc -2 ] is given by 

v, A r^sub /-R N i6o(161//m) i? W L(161/im) \~ 

* Lm=i «D — [-^— — — j , (1) 



-5 - 



where D is the distance to the SQ of 94 Mpc, A is the aperture area at the shocked region (260 
kpc 2 ), -Rni60 an d -Rwl are the relative responses of the N160 and WIDE-L bands at 161 /im, 
respectively, and A^niro and Az/wt , are the effective bandwidths of the N160 and WIDE-L, 
respectively. From lKawada et al. I (120071 ). -Rni6ch Rwl, Az/ N160 , and Az/ WL are taken as 0.96, 
0.59, 0.4 THz, and 0.8 THz, respectively. Thus, £l [cii] is estimated to be (1.0l ; 5 ) x 10 39 er § 
sec -1 kpc -2 , wh ich is comparable t o the H2 line luminosity surface density £l H2 of 2 x 10 39 
erg sec -1 kpc~ 2 (jCluver et al. Il2010l ). By using the dust emission at 160 /im, the EW of the 
[C II] line is estimated to be ~10 /jm. The EW is about 10 times larger than that in n earby 



galaxies (the mean of 0.7 /im with the dispersion of 0.2-3.0 /im, iBoselli et al. I (120021 ) ). We 



calculated the contribution of the [C II] line emission to the fluxes in the WIDE-L and N160 
bands by considering the constant dust emission in units of F v over the two bands. We 
obtained 8 % to WIDE-L band and 23 % to N160 band. 



3.2. Properties of far-IR dust emission 



To estimate the far-IR luminosity Lfir in the shocked region, the best-fit modified 
blackbody (the dash dotted line in Fig[3](a)) was integrated between 3 and 3000 /im to yield 
-Lfir of (5 ± 2) x 10 42 erg sec -1 . By dividing Lfir by the aperture area A, the luminosity 
surface density X_l fir is estimated to be (2.0 ± 0.6) x 10 40 erg sec -1 kpc -2 , which is in 
agreement wi th the far-IR lu minosity surface density (1.6 x 10 40 erg sec -1 kpc -2 ) of ISO 
observations (IXu et al. I 20031). In the shocked region, we calculate the dust mass according 
to Eq. (4) of lHildebrandl (119831 ) with the emissivity power-law index of unity. The dust mass 
surface density E Md is given by 



S A/d = 1.1 x 10 4 



2.0 x 10 40 erg sec -1 kpc 



-2 



22 K 



M kpc" 



(2) 



The dust temperature is set to be 22 K, which is derived from the SED fitting in Fig[3^a). 
From Eq. (jSJ), is estimated to be (1.1 ± 0.3) x 10 4 M Q kpc -2 . If we consider the 



Cluver et al. 



contri bution of the [C II] line emission to -F N160 , T, Md is reduced by ~20 %. 
( 120101 ) estimated the mass surface density of warm H 2 gas in the shocked region as (10 ± 
2) x 10 5 M & kpc -2 . CO observations of the shocked region suggest that the mass ratio of 
the cold H2 (< 50 K) to wa r m H2 gas masses is 2 w ith the CO-to-H2 conversion factor of 
our Galaxy ( jGuillard 1 120 101 ) . iGuillard et al. I (120091 ) estimated an upper limit of the mass 
surface density of H I gas in the shocked region as 5 x 10 5 M Q kpc -2 . Thus, the gas-to-dus t 
mass ratio is 170-210, which is in agreement with that of our Galaxy (jSodroski et al. Ill997l ). 
The estimated properties of the dust in the SQ are summarized in Table [TJ 
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Discussion 



4.1. Shock-powered [C 111158/um line emission from the shocked region 

To investigate the possibility of the luminous [C II] line emission from the shocked region, 
the C + abundance per hydrogen atom X c + is estimated with assumptions that the [C II] 
line emission comes from the warm H2 gas and is optically thin. Within these assumptions, 
Xq+ is given as 



X, 



c+ 



1.1 x 1(T 4 



J [CII] 



L kpc 



Sm h , M kpc 




2exp(-91 K/T) + n^/nn 
2exp(-91 K/T) 



(3) 



where T is the kinetic temperature of H2 gas, n^ 1 * is the critical density of the [C II] line 
emission for H I collisions, and nn is the nu mber density of H I At the shocked region, T and 
wh are 158 K and > 10 3 cm -3 , respectively (ICluver et al. 1120101 ) . is 3.0 x 10 3 cm -3 at 158 
K (jLanger et al. Il2010i ). Thus, Xq+ is estimated to be ~ 1 x 10~ 4 . Assuming that the carbon 
in the [CII] line emitting region is in singly ionized form, the carbon abundanc e in the shocked 



region is in agreement with that in an interstellar gas-phase (1.4 x 10 -4 , ICardelli et al. 
(119961 ) ) . Therefore, the luminous [C II] line emission from the shocked region is physically 
plausible provided that C + is the main carbon form in the warm H 2 gas. 



4.2. Origin of the cold dust emission from the shocked region 



The clear spatial correlations among the cold dust, H 2 , and X-ray emissions suggest two 
possibilities: radiative heating of dust grains in a low interstellar radiation field (ISRF) and 
collisional heati ng of dust grains in a hot plasma. To explain the formation of H2 gas in the 
shocked region, iGuillard et al. I (120091 ) introduced a model of the shock in an inhomogeneous 
gas phase. In their model, there are two critical assumptions to explain the formation of 
H2 gas in the collision age: dust survival in the shock and the Galactic gas-to-dust mass 
ratio. The spatial correlation between the maps of the cold dust emission and H2 rotational 
line strongly supports H 2 formation on dust grains. In addition, the observed gas-to-dust 
mass ratio in the shocked region is in agreement with the Galactic one. These facts indicate 
that the H 2 formation time scale should be shorter than the collision age. In the case of the 



multip hase medium, the cold dust emission comes from molecular gas clouds. IGuillard et al. 



(120101 ) modeled SEDs of dust emission associated with diffuse (the hydrogen column density 



Nft ~ 2 x 10 20 cm -2 ) or clumpy molecular gas (Nr ~7x 10 21 cm -2 ) in radiation fields. As 
discussed in Section 4.1, the [C II] line emission can come fro m molecular gas c l ouds. The 
optical depth of the [C II] line t\co] is expressed in Eq. (A4) of ICrawford et al. I (119851 ). By 
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using X c + of 1 x 10~ 4 and the velocity dispersion of 870 km s -1 ( ICluver et al. Il2010l ). Ah 
is estimated to be 2 x 10 24 cm -2 when rrciii becomes unity. Thus, the [C II] line emission is 
optically thin for both cloud models. To estimate color temperatures of col d dust for both 



cloud models, flux densities at 90 and 160 /j,m are obtained from Fig. 7 in iGuillard et al. 



(120091 ). Assuming the modified blackbody spectrum with the emissivity power-law index of 
unity, the color temperature is estimated to be 24 K for the diffuse molecular gas model and 
20 K for the clumpy one. We investigate the spatial distribution of the color temperature 
in the SQ (see Fig. 0J by using the flux densities at 90 and 160 //m that are above the 
3a sky background noise. Errors in the color temperature are estimated as 2.0-2.5 K. The 
color temperature in the shocked region is about 20 K, which favors the temperature for 
the clumpy molecular gas model rather than the diffuse one and significantly lower than 
that in the surrounding galaxies. Thus, it is possible that cold dust in the shocked region is 
radiatively heated by the ISRF. 

The spatial correlation between the maps of the cold dust emission and X-ray suggests 
an alternative possibility for the origin of cold dust emission in the shocked region. If dust 
grains are exposed to the hot plasma for the collision age (~ 5 x 10 6 yr), grains smalle r 
than 0.1 fim in radius must have been destroyed in ~ 5 x 10 6 yr (IGuillard et al. II2009I ). 
We calculated the temperatures of dust grains with radii a > 0.1 /xm h eated collisionally 



by ambient plasma electrons, based on the model of iDwek I (119861 . 119871 ). We assumed the 
size distribution dn/da oc a~ a with a = 2.5-3.5 over the range 0.1 < a < 1 fim for either 
silicate or graphite grains. The density and temp erature of the electron s are taken to be 
n e = 0.02 cm -3 and T e = 6 x 10 6 K, respectively ( ITrinchieri et al. 1120031 ). As a result, we 
derive the dust temperature of 20-22 K, which is insensitive to a specific choice of a or grain 
composition. The derived temperature is in a range similar to that observed in the shocked 
region. Thus we cannot distinguish between radiative heating and collisional heating of dust 
grains solely from the dust temperature. 

However, the similarity in the spatial distribution between the far-IR (cold dust and 
[C II] line) emission and H2 emission s trongly suggests tha t H2 molecules form on dust 
grains and coexists with C + . Moreover, iNatale et al. I (2010J ) showed that the collisionally 
heated dust emission should be minor contribution to explain the X-ray temperature map. 
Therefore, it is likely that the far-IR dust emission arises mostly from radiative heating of 
cold dust in clumpy molecular gas. These facts support the multiphase medium model for 
the SQ. 
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5. Summary 

We observed the SQ with the FIS aboard AKARI in four far-IR bands. The N160 
band image shows single peak emission in addition to the structure extending in the North- 
South direction along the shock ridge as seen in the WIDE-L band, H 2 emission and X-ray 
emission. Whereas most of the far-IR emission in the shocked region comes from the cold 
dust component, the [C II] line emission whose luminosity is comparable to that of the warm 
H 2 gas can significantly contribute to the single peak emission in the N160 band. It can be 
explained that the [C II] line emission comes from the warm H2 gas. In the shocked region, 
the observed gas-to-dust mass ratio is in agreement with the Galactic one. These indicate 
that H 2 molecules are produced on the grain surface and their formation timescale should 
be shorter than the collision age. The color temperature in the shocked region (~ 20 K) 
is significantly lower than that in the surrounding galaxies (~ 30 K). It is considered that 
the cold dust emission comes mostly from radiative heating in dense and clumpy gas clouds. 
Our results support the scenario of the multiphase medium of the pre-shocked gas in the 
IGM. 

We are grateful to the referee for providing us very important comments and correc- 
tions on our calculation. The present work is based on observations with AKARI, a JAXA 
project with the participation of ESA. The Digitized Sky Surveys were produced at the 
Space Telescope Science Institute under U.S. Government grant NAG W-2166. The images 
of these surveys are based on photographic data obtained using the Oschin Schmidt Tele- 
scope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed into 
the present compressed digital form with the permission of these institutions. This work was 
supported in part by Grant-in-Aid for Young Scientists by MEXT (21740139). 
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Table 1: Properties of the far-IR dust emission at the shocked region 





^FIR 




Gas-to-dust 


(K) 


(erg sec -1 kpc~ 2 ) 


(M kpc" 2 ) 


mass ratio 


22±1 


(2.0 ±0.6) x 10 40 


(1.1 ±0.3) x 10 4 


170-210 




Fig. 1. — Composite image of the Stephan's Quintet (HCG92): AKARI 7 /j,m (red), Optical 
(green) , and X-ray (blue) . Optical and X-ray images are retrieved from Digitized Sky Survey 
and XMM-Newton Science Archive (observation ID: 0021140201, observation date: 2001-12- 
07), respectively. The SQ is composed of five galaxies (NGC 7317, NGC 7318a NGC 7318b, 
NGC 7319, and NGC 7320c). NGC 7320 is an unrelated foreground object jMoles et al. 
19971 ). The box shows the boundary of the maps in Fig. [21 
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(f) Contours: X-ray 
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Fig. 2.— Four-band images of the SQ in the N60 (top-left), WIDE-S (top-middle), WIDE-L 
(bottom- left), and N160 (bottom- middle) bands. The center wavelengths of the four bands 
are 65 /im for N60, 90 /im for WIDE-S, 140 /im for WIDE-L, and 160 /im for N160. The 
contours are linearly spaced from 10 % to 99 % of the peak brightness with a step of 3 %. 
The peak brightness that includes the sky background is 22 MJy sr -1 (N60), 17 MJy sr _1 
(WIDE-S), 22 MJy sr" 1 (WIDE-L), and 27 MJy sr" 1 (N160). In each image, the PSF size 
in FWHM is shown in the lower left corner. The upper-right and bottom-right panels show 
the images of WIDE-S and N160 bands overlaid on AKARI 7 /im and XMM/Newton X-ray 
contours, respectively. The three apertures for photometry of the shocked region (r = 20"), 
NGC 7319 (r = 30"), and NGC 7320 (r = 30") are shown in the panels of (c) and (f). As a 
reference, the shock ridge is shown by the solid white line in the panel (f). 
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Fig. 3.— (a) SEDs at the shocked region, NGC 7319, and NGC 7320. The positions of the 
apertures are shown in Figs [2(c) and (f). The dashed, dotted, and dash-dotted lines show 
the single temperature modified blackbody spectrum with the emissivity power-law index of 
unity for NGC 7319, NGC 7320, and the shocked region, respectively, (b) The N160 image 
shown after subtraction of the WIDE-L image. The contours show XMM/Newton X-ray 
contours that are the same as those in Fig. [2(f). 
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Fig. 4. — Color temperature map of the SQ. The contours show XMM/Newton X-ray con- 
tours that are the same as those in Fig. MS)- 



